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Abstract. We present RXTE and BeppoSAX observations of the Intermediate Polar V709 Cas acquired in 
1997 and 1998 respectively. The X-ray emission from 0.1 to 30keV is dominated by the strong pulsation at the 
rotational period of the white dwarf (312.8 s) with no sign of orbital or sideband periodicity, thus confirming 
previous ROSAT results. However, we detect changes in the power spectra between the two epochs. While the 
second harmonic of the spin period is present during both observations, the first harmonic is absent in 1997. An 
increase in the amplitude of the spin pulsation is found between 1997 and 1998 together with a decrease in the 
X-ray flux. The average X-ray spectrum from 0.1 to 100 keV is well described by an isothermal plasma at ~ 
27keV plus complex absorption and an iron Kc fluorescent line, due to reflection from the white dwarf surface. 
The rotational pulsation is compatible with complex absorption dominating the low energy range, while the high 
energy spin modulation can be attributed to tall shocks above the accreting poles. The RXTE spectrum in 1997 
also shows the presence of an absorption edge from ionized iron likely located in the pre-shock accretion flow. 
The variations along the spin period of the partial covering absorber and of reflection are compatible with the 
classical accretion curtain scenario. The variations in the spin pulse characteristics and X-ray flux indicate that 
V 709 Cas experiences changes in the mass accretion rate on timescales from months to years. 
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1. Introduction 

Intermediate Polars (IPs) are a subclass of magnetic 
Cataclysmic Variables (mCVs) consisting of a weakly 
magnetized (up to a few MG) and asynchronously rotat- 
ing (Prot < Porb) whitc dwarf (WD) which accretes ma- 
terial from a late type, main sequence, Roche-lobe filling 
secondary star. 

The magnetic field influences the details of the flow down 
to the WD poles. In low field IPs, an accretion disc, trun- 
cated at the magnetospheric radius, may form and the ma- 
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terial flows from the inner edge of the disc along the fleld 
lines in an arc-shaped curtain towards the polar regions 
of the WD (Rosen et al. 1988). A disc- less (or stream- fed) 
accretion is expected in systems with higher fields, where 
material accretes directly from the stream flowing from the 
secondary star. A mixture of stream-fed and disc accre- 
tion are also possible, where matter from the stream over- 
passes the disc (disc-overflow, Hellier et al. 1995). Different 
configurations (disc-fed and disc-overflow) can be present 
at different epochs in the same system, due to long term 
changes in the mass accretion rate (de Martino et al. 1995; 
Buckley et al. 1996; de Martino et al. 1999). 
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The inflow of accreting material onto the polar caps pro- 
duces a stand-ofi^ shock, below which material cools via 
thermal Bremsstrahlung (kT ~ 5 — 30 keV) . Hard X-rays 
are highly absorbed by cold matter (column densities up 
to 10^"^ cm~^) within the accretion flow (Ishida et al. 1994) 
and are expected to be reflected by the WD atmosphere 
(Beardmore et al. 1995; Done et al. 1995; Matt 1999). 
Because of their asynchronous rotation, IPs show a wide 
range of periodicities at the WD spin (oj), the orbital 
(n) and sideband frequencies (Warner 1986; Patterson 
1994, Warner 1995), whose amplitudes can be wavelength 
dependent (de Martino 1993). The X-ray periodicities 
provide direct information on the actual accretion mode 
(Wynn & King 1992; Norton et al. 1996), disc-fed systems 
being dominated by the WD rotational period. The addi- 
tional presence of periodicities, at the beat (or synodic) 
u> — n and orbital (17) frequencies are an indication of 
disk-overflow accretion. A pure stream-fed system should 
show pulses at the synodic frequency only (Hellier 1999). 
To date only V2400 Oph is a confirmed stream-fed IP 
(Buckley et al. 1997). Furthermore disc-fed systems show 
different rotational pulse characteristics: double-peaked 
X-ray pulses, generally observed in fast rotators such as 
YYDra, can originate from the emission of two accret- 
ing poles with tall shocks or large footprints (Norton et 
al. 1999), whilst single-peaked curves, observed in most 
slow rotators such as FO Aqr, can instead be produced 
by small accreting regions. It has been suggested (Norton 
et al. 1999) that such differences are due to the effect of 
the magnetic field strength, the rapid rotators containing 
low field WDs (Norton et al. 1999). Unfortunately, for IPs 
such a proposal is difficult to confim observationally, due 
to the lack of detectable optical-IR polarization in most 
cases and of spectral features of the WD. An outstand- 
ing exception is V709 Cas (Bonnet-Bidaud et al. 2001), 
the first IP known to date to display the WD atmospheric 
features, from which a magnetic field strength lower than 
3 MG has been estimated. 

V709 Cas (RX J0028.8+5917) was identified as an IP 
from the ROSAT AU Sky Survey (Haberl & Motch 1995; 
Motch et al. 1996). The X-ray emission was found to be 
hard (HR2=[1. 0-2.4] keV - [0.4-1.0] keV/[0.4 - 2.4] keV) 
= -1-0.21) and pulsed at a period of 312.8 s in pointed 
follow-up PSPC and HRI ROSAT observations (Haberl & 
Motch 1995; Norton et al. 1999) (see Table 1 for the history 
of pointed X-ray observations). Searches in the HEAO- 
1, UHURU and Ariel V catalogues lead to the identifica- 
tion of V 709 Cas with the previously detected hard X- 
ray sources IH 0025+588, 4U 0027+59 and 3A 0026+593 
(Motch et al. 1996). The X-ray timing analysis shows the 
presence of the dominant 312.8s pulsation and its first 
and second harmonics, with no sign of orbital and side- 
band variabilities (Norton et al. 1999). Two possible or- 
bital periods at Po=5.4h and Pf2=4.5h have been in- 
ferred from low temporal resolution optical spectroscopy 
(Motch et al. 1996), each being the one day alias of the 
other. The former has been recently confirmed spectro- 
scopically (Bonnet-Bidaud et al. 2001) and photometri- 



cally (Kozhevnikov 2001). This suggests that V709 Cas is 
a classical IP: a hard X-ray emitting and disc-accreting 
mCV. 

In this paper we present X-ray observations of V 709 
Cas acquired with the RossiXTE (RXTE) and the 
BeppoSAX satellites. The latter also provides the first si- 
multaneous soft and hard X~ray data which allow a final 
characterization of its X-ray spectrum and its variability. 

2. Observations and data reduction 

2.1. The RXTE data 

V 709 Cas was observed with RXTE (Bradt et al. 1993) 
between March 26-28, 1997, over a 2.3 day interval. As the 
source was not significantly detected with the HEXTE in- 
strument, operating in the energy range 20-200 keV, we 
limited the analysis to the PCA instrument, consisting of 
an array of 5 proportional counter units (PCU) sensitive 
in the 2.5-60.0 keV band. The effective on-source expo- 
sure time, after standard screening, was 40 ks. V709 Cas 
was detected at a net count rate of 7.04+0.01 ctss~^ per 
PCU in the 2.5-25 keV band. We have estimated the PCA 
background using a faint source model appropriate for the 
epoch of the observation. The RXTE data are telemetred 
in two standard, and up to four additional, data modes. 
Two separate observation modes were used for spectral 
and timing analysis. For the production of spectra, we 
used Standard 2 mode data, which have a time bin size 
of 16 seconds. In order to improve the signal to noise, we 
used only data from the top layers of the PCUs. This re- 
duces the net signal by 20% and the estimated background 
by roughly 50%. Since spectral uncertainties, particularly 
above lOkeV, are dominated by systematic uncertainties 
in the background model, this improves the quality of the 
net spectrum. For timing analysis, we have used the Event 
mode data, with a 8 ms time resolution and 32 PHA chan- 
nels. As no selection by PCU or PCU layer is possible with 
this data mode, all layers from all PCUs were used. All 
analyzed data refer to the effective 2.9-24.2 keV band. 
The average background level in the five PCU units is 
77.1 ctss^^ from all layers, while it is 40.3 ctss^^ from the 
top layer. 

2.2. The BeppoSAX data 

A later observation with the BeppoSAX satellite (Boella 
et al. 1997) was carried out between 1998 July 5-7 
with the set of the co-aligned Narrow Fields Instruments 
(NFI), used as prime pointing instruments, covering 
the wide energy range 0.1 - 300 keV. The source 
was detected with the LECS (0.1-lOkeV), the MECS 
(1.3-lOkeV) and PDS (15-300keV) detectors at net 
count rates of 0.243+0. 003ctss-i, 0.497+0. 002ctss-i 
and 0. 346+0. 028ctss~^ respectively. Effective on-source 
exposure times have been 35.3 ks for the LECS, 86.9 ks 
for the MECS and 79.5 ks for the PDS detector (two units, 
see below). 
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Table 1. History of pointed X-ray observations of V709 Cas. 



Date 



Instrument 



Energy Band 
(keV) 



Count Rate 

ctss""*" 



Flux 

-2 



erg cm 



1992, Jul. ROSATPSPC 0.1-2.4 

1997, Mar. RXTEPCA 2-10 

1998, Feb. ROSATHRI 0.1-2.4 
1998, Jul. BeppoSAX MFCS 2-10 



0.60 
5.20 
0.26 
0.50 



8.8x10"^^ 
6.6x10"" 
~ 1x10"" 
4.5x10"" 



Source 



a: Haberl & Motch (1995). 
b: This work. 

c: Norton ct al. (1999). Note: The X-ray flux is estimated using the count rate conversion facility PIMMS-HEASARC aX 

cm 



GSFC-NASA, adopting kT=10keV and Nh = 9.8 x 10^" 



The LEGS light curves and spectra have been extracted 
using a circular region with a radius of 8 arcmin centred 
on the source, whilst MEGS light curves and spectra have 
been extracted with a smaller 4 arcmin radius. While the 
whole LEGS band will be used for timing analysis, the 
spectra from this instrument have been analyzed only be- 
low 4 keV due to calibration problems at higher energies. 
For both intruments, the background has been measured 
and subtracted using the same detector regions during 
blank field pointings. The measured background in the 
LEGS (2.74xl0"2ctss-i) and MEGS (6.95x10--'' ctss"!) 
is much lower than the source count rate. 
The PDS is a coUimated instrument which monitors the 
background continuously switching two (of the four) de- 
tectors with a dwell time of 96 s. The PDS light curves and 
spectra have been extracted using a standard routine pro- 
vided by the BeppoSAX Data Genter (SDG). We have con- 
servatively extracted the light curves in the energy range 
of 15-70 keV, above which the source is hardly detected. 
The average net count rate in this range is 0.270± 0.020 
ctss~^. The procedure allows rejection of particle back- 
ground events, as well as spikes caused by single particle 
hits, which produce fluorescent cascades which are mainly 
recorded below 30 keV. The background spectrum is eval- 
uated for each of the two half arrays accumulating the 
off-source spectra. The background light curves are con- 
structed by linear interpolation between two off-source 
pointings in each array. Then the spectra and light curves 
from the two arrays are merged together to construct the 
background spectrum and light curve for final subtraction. 
The average background count rate is 10.31 ctss"^. We 
also performed alternative ways to subtract background, 
accounting for collimator dwell time effects in order to 
construct the net PDS light curve for periodicity search 
purposes. We detected no significant modulation in the 
15-70 keV energy range (see sect. 3). 

The history of the pointed X-ray observations of V 709 
Gas is reported in Table 1. 



3. Timing analysis 

A search for periodicities has been performed on light 
curves extracted from both data sets in similar energy 
bands. The range 2.9-9.8 keV has been used for the RXTE 
PGA data, and the 1.3-lOkeV band for the BeppoSAX 
MEGS data. In Fig. 1. the RXTE and BeppoSAX hght 
curves with a temporal resolution of 32 s are shown. A 
short term periodicity is clearly apparent, with no indi- 
cation of variability (periodic or not) on a longer (hours) 
timescale. 

Fourier analysis has been performed on PGA 4 s and 
MEGS 5 s binned light curves using the DFT algorithm 
(Deeming 1975). The corresponding power spectra, shown 
in Fig. 2 (lower panels), reveal a strong peak at the known 
312.8 s period at both epochs. While the first and sec- 
ond harmonics are observed in the MEGS spectrum, the 
PGA spectrum lacks power at the first harmonic. In or- 
der to remove the windowing effects of unevenly sampled 
data (see spectral windows inserted in the lower panels of 
Fig. 2), the GLEAN algorithm (Roberts et al. 1987) has 
been used, adopting a gain of 0.1 and 500 iterations. From 
the GLEANed power spectra, shown in Fig. 2 (upper pan- 
els), neither low frequency (orbital) periodicity nor the 
beat or other orbital sidebands can be detected. 

From the MEGS and PGA GLEANed spectra the ampli- 
tude of the main signal is A^ = 0.108 ± 0.001 ctss^^ and 
Aj^ = 0.82 ± O.Olctss"^ respectively, corresponding to a 
fractional amplitude of 23% and 18% at the two epochs. A 
larger amplitude in the July 1998 BeppoSAX data is also 
found for the second harmonic, Aa^ ~ 0.034±0.002 cts s~^ 
(fractional amplitude 7.5%), with respect to the March 
1997 data, A^^ = 0.33 ± 0.01 cts s"^ (fractional ampli- 
tude of 6.5%). This also translates in the lack of detec- 
tion of the first harmonic in the RXTE data while in the 
BeppoSAX observations the first harmonic is definitively 
detected, with an amplitude of A2cj — 0.019±0.002ctss-^ 
(fractional amplitude 3.8%). This fullfils the 4ct criterium, 
with a — 4.3 x 10~''ctss~^ being the average amplitude 
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Fig. 1. The net RXTE PCA (2.9-9.8 keV) (left panel) and BeppoSAX MECS (1.3-10 keV) (right panel) 32 s binned 
light curves of V 709 Cas in March 1997 and July 1998 respectively. 



from the residual CLEANed spectrum. Here we note that 
the ROSAT HRI observation carried out in February 1998 
(Norton et al. 1999), and hence about four months earlier 
than the BeppoSAX pointing, also shows the presence of 
the first and second harmonics. However, the ratios of the 
power at the fundamental frequency to the first and to the 
second harmonics in the BeppoSAX LEGS (0. 1-2.4 keV) 
range are about 90 and 10 respectively , a factor of three 
and two larger than those derived from the ROSAT HRI 
observations. 

A sinusoidal fit to the MECS light curve using three 
frequencies corresponding to the fundamental, first and 
second harmonics yields to a period of 312.77 ± 0.01s. 
Instead, for the PCA light curve two sinusoids, cor- 
responding to the fundamental and second harmonic, 
have been used, which give P^ — 312.746 ± 0.003 s. 
The present determinations are compatible within er- 
rors with the previously determined 312.78±0.03s period 
from the ROSAT data (Haberl & Motch 1995; Norton 
et al. 1999). Times of maxima for each data set are 
derived: HJO^ax^ 2450533.647454(3) and HJDmax^ 
2451000.32346(2). The accuracy of the derived ephemeris, 
however, does not allow to bridge the gap between the 
RXTE and BeppoSAX data sets. 

The shape of rotational pulse has been inspected in dif- 
ferent energy bands (Fig. 3). The RXTE data in the 2.9- 
4.7 keV, 4.7-9.8 keV and 9.8-24.2 keV bands (Fig. 3, left 
panel) have been folded in 50 bins along the rotational 
period using the corresponding time of maximum. The 
BeppoSAX data instead have been extracted in the 0.1- 
2.4 keV LECS range (for a direct comparison with the 



ROSAT hght curve see Fig.4), in the 2-5 keV and 5-lOkeV 
MECS ranges as well as in the 15-70 keV PDS band. The 
spin light curves, also folded in 50 bins (except for the PDS 
where a larger bin size (21 s) has been used), are shown in 
Fig. 3 (right panel). At both epochs, the rotational pulse 
is not sinusoidal and shows a broad structured maximum 
(extending 0.46 in phase). The spin pulse as observed by 
RXTE has relatively symmetrical rise and decay. On the 
other hand, the BeppoSAX light curve is highly asymmet- 
rical with a slow rise and a steeper decay. The maximum 
is not flat with count rates increasing until = 0.12. A dip 
feature with fractional depth ^ 8% can be recognized in 
the hard MECS bands at 0=0.9. This feature does not oc- 
cur at mid-maximum and it is less pronounced in the soft 
band. This feature could also be present in the RXTE light 
curve, but its depth is only marginally significant. The 
changes in the shape of the spin modulation between the 
two epochs then explain the presence of the first harmonic 
in the BeppoSAX data and the lack of such component in 
the RXTE data. 

The spin pulsation decreases in amplitude from low to high 
energies: the pulsed fraction in the RXTE data changes 
from 46% in the 2.9-4.7keV range to 32% in the 9.8- 
24.2 keV band, and in the BeppoSAX data it decreases 
from 86% in the 0. 1-2.4 keV, to 56% in the 2-5 keV and to 
52% in the 5-10 keV band. In the PDS range, the variabil- 
ity is not statistically significant. The light curves reveal 
an energy dependence of the pulse shape changing from 
a strongly asymmetrical (at low energies) to a broad (at 
high energies) structured shape. A comparison of the low 
energy pulsation in the 0.1-2.4keV (ROSAT-like) band 
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Fig. 2. The PCA (left panel) and the MECS (right panel) power spectra of V709 Cas. Lower panels: The DFTs 
show a strong signal at 276.26 d^^. In July 1998 both first and second harmonics are present, while in March 1997 the 
first harmonic is not detected. The spectral windows are shown in the inserted panels. Upper panels: The removal of 
windowing effects is shown in the CLEANed spectra. 
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Fig. 4. The ROSAT HRI rotational pulsation in the 0.1- 
2.4 keV band, observed by Norton et al. (1999) has been 
folded at the period of 312.77 s adopting a bin size of 6.26 s. 
Phase zero is arbitrary (see text). 



observed in July 1998 with the February 1998 ROSAT 
HRI rotational curve (Norton et al. 1999) shows a differ- 
ent morphology. The ROSAT spin pulsations is reported 
in Fig. 4 with an arbitrary phasing since the BeppoSAX 
ephemeris is not accurate enough. The fractional ampli- 
tude of the ROSAT spin pulse is about half that ob- 
served by BeppoSAX. Also, the broad maximum in the 
BeppoSAX data still covers the same phase range but the 
"double-peaked" shape seen in the ROSAT data, with a 



secondary minimum centred on the broad maximum, has 
essentially disappeared. The dip in the high energy light 
curve does not seem to have a counterpart in the soft X- 
ray pulse. This feature is not centered on the maximum 
as observed in the ROSAT light curve. Although a direct 
comparison of the RXTE and ROSAT light curves should 
not be performed due to the different energy coverages, the 
flat-topped RXTE light curve resembles more closely the 
ROSAT spin pulsation with two almost identical maxima. 
Again, the detection of an increasing power of the first 
harmonic in the RXTE, ROSAT and BeppoSAX observa- 
tions accounts for such differences. 

The energy dependence of the rotational pulse manifests 
itself in the hardness ratios in both soft and hard X-ray 
bands. In Fig. 5 (right panel), the LEGS (2-5keV/0.1- 
2keV) and MECS (5-10keV/2-5keV) ratios show a hard- 
ening at rotational minimum and a softening at rotational 
maximum. This behaviour is also observed in the RXTE 
PCA data (Fig. 5, left panel), which also shows a negligible 
energy dependence of spin pulses above lOkeV. A hard- 
ening at pulsation minimum is a typical behaviour in IPs, 
which is produced by the larger photoelectric absorption 
when viewing along the accretion curtain. We note that 
the BeppoSAX hardness ratios do not reveal changes dur- 
ing the maximum phases, which would be expected if the 
double-humped maximum is due to enhanced absorption 
effects within one accretion column. As it will be discussed 
later, the presence of two accreting poles can account for 
this behaviour. At higher energies a different mechanism 
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Fig. 3. Left panel: The RXTE light curves in the 2.9-4.7keV, 4.7-9.8 keV and 9.8-24.2 keV ranges (from bottom to 
top) folded along the rotational period using the ephenieris reported in the text. Right panel: The BeppoSAX net spin 
light curves in the energy ranges (from bottom to top): 0. 1-2.4 keV, 2-5 keV, 5-10 keV and 15-70 keV. A 6.26s bin size 
has been used except for the PDS (15-70 keV) light curve for which a 21 s bin has been adopted. 



responsible for the spin pulsation could be at work as it 
will be discussed in sect. 4 and 5. 



4. Spectral Analysis 

An analysis of the time-averaged X-ray spectra of both 
RXTE and BeppoSAX observations has been performed 
in order to estimate the basic parameters of the emit- 
ting region, in particular the temperature of the optically 
thin plasma, interstellar and possible circumstellar ab- 
sorptions and to obtain information on the reprocessing of 
the hard X-rays from the WD surface. The spectral anal- 
ysis has been performed with the XSPEC fitting package. 
Henceforth, all quoted errors on spectral parameters refer 
to 90% confidence level for the parameter of interest. The 
BeppoSAX spectra have been fitted between 0.3-70 keV, 
while the RXTE ones between 3.6-20 keV. 



4.1. The BeppoSAX spectra 

A spectral fit consisting of an optically thin plasma 
(MEKAL model (Mewe et al. 1995) in XSPEC) with a 
single component absorber does not produce a satisfac- 
tory fit (x^/dof=243/114), the temperature being as high 
as 60keV. The hydrogen column density is 2.6 x 10^^ cm^^ 
which is lower than the total galactic hydrogen column 
density (4.5x10^^ cm^^) in the direction of the source. 
The residuals show an excess of counts in the low energy 
portion of the spectrum, likely due to a more complex 



absorption, as well as around ~6 keV. Hence we have in- 
cluded a Gaussian Kq, irone line centred at 6.4 keV and 
an additional (partial) absorber. The presence of this ab- 
sorber is also suggested by the energy dependence of the 
spin pulse and the strong (see below) iron line. We then ob- 
tain a much better fit {xl — 1.17) (see Table 2, model 1). A 
lower temperature is derived, namely kT = 42 keV which 
reflects the effects of a complex model fltting. Here we also 
note that a multi-temperature plasma model (CEMEKL 
in XSPEC, which is built from the MEKAL code, adopt- 
ing a power-law temperature profile (Singh et al. 1996)) 
gives a worse fit {xl =1-5), especially for the ionized iron 
lines and the low energy portion of the spectrum. The hy- 
drogen column density, Nh = 8.7 x 10^° cm^^, is close to 
that derived from ROSAT data (Nr = 9.8 x 10^° cm'^ 
for a fixed temperature of kT=10keV, Haberl & Motch 
1995). This can be regarded as an upper limit to the in- 
terstellar column density. For the partial absorber we find 
Nh = 2.92 x 10^^ cm~^ and a covering fraction Cp = 0.29. 
The equivalent width of the 6.4 keV fluorescent line is 
E.W.=218eV. The metal abundance, left free to vary, 
is just consistent with the solar value. Given the high 
temperature of the single temperature plasma model, the 
metal abundance is only sensititive to the iron abundance. 
We have adopted a relative iron abundance by number of 
4.7x10"'"' (Anders & Grevesse 1989). The observed flux 
in the 2-lOkeV range is 4.5xl0~^^ergcm~^ s"^, a factor 
2.5 higher than that extrapolated from the ROSAT 1992 
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Fig. 5. The behaviour of hardness ratios in the soft and hard X-ray bands along the rotational period. Left panel: 
the RXTE PCA 5-10 keV and 3-5 keV ratios (lower panel) and the 10-24 keV and 5-10 keV ratios (upper panel). 
Right panel: the BeppoSAX MECS 5-10 keV and 2-5 keV ratios (upper panel) and the LEGS 2-5 keV and 0.1-2 keV 
flux ratios (lower panel). Both show a spectral hardening at rotational minimum. The effect decreases at increasing 
energies. 



PSPC data, adopting kT=10keV, and it is similar to other 
measurements by previous satellites (Motch et al. 1996). 

A strong K^ iron line would originate by fluorescence ei- 
ther from the surface of the WD or from the pre-shock 
accretion column, or both. To account for the large equiv- 
alent width, the column density of cold matter should be 
as high as 2 X lO^'^cm^^ (Inoue 1985), while a factor of 
~10 lower is found for the partial absorber. Hence it seems 
likely that the bulk of the fluorescent line comes from the 
WD surface. This suggests the presence of the contribu- 
tion of a Compton reflected continuum from the white 
dwarf which is expected to go along with the fluorescent 
iron hne (Matt et al. 1991; Done et al. 1992). Therefore, 
the introduction of a Compton reflection continuum is 
physically plausible, and turns out to be significant at 
the 99.7% confidence level, according to the F-test (see 
Table 2, model 2). Both temperature and fiuorescent line 
parameters change when adding such component: a lower 
temperature (27keV) and E.W.=179eV are derived from 
this fit. Hence, neglecting the reflected contribution the 
temperature and the strength of the fluorescent line are 
overestimated. The grand average combined LEGS, MECS 
and PDS spectra fltted with such a model are shown in 
Fig. 6. The relative normalization of the reflection com- 
ponent, represents the solid angle subtended by the cold 
matter in units of 2tt steradians for an average value of 
cos 9 — 0.5 where 9 is the angle between the line of sight 
and the normal to the reflecting surface. The value 0.92 is 



compatible with refiection of photons from the WD sur- 
face which subtends an angle of about 27r. The iron and 
metal abundances (linked together in this fit) are within 
errors consistent with solar values. Furthermore the par- 
tial absorber covering fraction and column density are, as 
expected, the same as in the previous fit without reflec- 
tion. 

The variability in the hardness ratios is a strong evidence 
that spectral changes occur at the rotational period. A 
spin resolved spectral analysis has been performed using 
only the LECS and MECS data, due to the lack of infor- 
mation on variability from PDS data. Spectra have been 
extracted in 10 rotational phase bins. In order to estimate 
changes in the spectral parameters, a fit with a composite 
model consisting of a MEKAL emission plus the inter- 
stellar column density, a partial absorber and Gaussian 
line has been performed. With the exclusion of the PDS 
data the temperature of the optically thin gas has been 
fixed at 30keV, which is the value obtained from fitting 
the average spectrum without the PDS data to model 1. 
We note that the fit quality does not substantially change 
if kT=40 keV is assumed. As a first attempt both cover- 
ing fraction Cp and Njj of the partial absorber and the 
Gaussian line normalization have been left free to vary, 
with the remaining parameters fixed at the values found 
for the phase-averaged spectrum (model 1). In these fits, 
the partial covering parameters are less well defined, es- 
pecially Nh, which is essentially constant within errors. 
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Fig. 6. Left Panel: The grand average BeppoSAX LEGS, MEGS and PDS spectra fitted with an optically thin plasma 
at kT=27keV together with an interstellar and partial covering absorbers and a Gompton reflected continuum and 
Gaussian centred at 6.4 keV (see text). Right Panel: The X-ray spectrum corresponding to the composite model is 
displayed in physical units. 

Table 2. Best fit parameters for the phase-averaged X-ray spectra of V 709 Gas 



Data 


# 


Nh 
(lO^o cm-2) 


Nh" 
(10^2 cm-2) 


c^ 


kT" 
(keV) 


Az' 


E.W.<= 

(eV) 


R-rcf 


T 


Eo 
(keV) 


xVd.o.f. ixt) 


SAX 
SAX 


1 
2 


8.7l:li 

8.011° 


9 q+1.6 

^•^-0.9 
o O+0.9 
^"^-0.7 


0.29l«:« 
0.301H? 


42+? 
27tl 


1.29l°:i 

r, yo + 0.28 
'-'•''^-0.20 


218t- 
179tll 


0-9tl;t 






131/111 (1.18) 
121/110 (1.10) 


RXTE 

RXTE 


3 

4 




1.61'^:? 

l-7tV2 




36+1 
26tl 


n co+0.09 

0.40_o.i7 


216tll 
198+i« 


3+"-2 


'-'•J^J-0.02 

n 11 +0-02 


8.0t°i 


45/37 (1.23) 
40/36 (1.11) 



" Column density of the partial absorber. 
*" Covering fraction of the partial absorber. 
■^ Plasma temperature. 

Metal abundance in units of the cosmic value (Anders & Grevesse 1989). 
* Equivalent width of the 6.4 keV fluorescent iron line. 
^ Relative normalization of the reflection component (see text). 



Hence fits have been performed again varying only Gp 
and the Gaussian normalization, with Njj of the partial 
absorber fixed at the value found for modell. The results 
are shown in Fig. 7, where the changes in the covering frac- 
tion (a factor ^ 1.6) indicate an increase in absorption at 
rotational minimum. The E.W. of the fluorescent line also 
increases at rotational minimum by a factor ^ 2.3. 

An attempt to search for variations in the reflection 
component has been performed using the spectra at 
rotational maximum and minimum, keeping constant 



the temperature, abundances, interstellar and partial 
covering hydrogen column densities at the values found 
from flt with model 2, in Table 2. The variations of the 
relative normalization Rrcfi in these flts are consistent 
with the iron line E.W. behaviour shown in Fig. 7, being 
lower at rotational maximum (Rrcfi < 0.2) and larger at 
minimum (Rrofi = 0.7 ± 0.4). 
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Fig. 8. Left Panel: The grand average RXTE spectrum fitted with an optically thin plasma at kT=26 keV together 
with an absorber with Nh — 1-7 x 10^^ cm^^, an absorption edge with r — 0.11 and Eq — 8.0 keV, and a Compton 
reflection continuum plus a Gaussian line centred at 6.4 keV (model 4). Right Panel: The X-ray spectrum corresponding 
to the composite model is displayed in physical units. 



4.2. The RXTE spectra 

The time-averaged RXTE spectrum was fitted first with a 
MEKAL model, together with a Gaussian line centred at 
6.4 keV, and two absorbers to account for the interstellar 
and partial covering but the fit quality is bad (xt = 5.97). 
The residuals show a strong excess of counts around 5 keV 
and a depression around 8 keV. The presence of this fea- 
ture has been checked against instrumental (background 
subtraction) effects, and in the individual spectra of all five 
PCU units. We interpret this feature as a signature of an 
absorption edge. Furthermore the interstellar column den- 
sity results are negligible, as expected by the limited ex- 
tension to low energies. The covering fraction of the partial 
absorber results instead close to unity (total absorber). 
Hence the fits have been performed again neglecting the 
interstellar absorption and introducing an additional com- 
ponent (e~'^(^/^°) at E > Eo) in order to account for 
the absorption edge. As shown in Table 2 (model 3), the fit 
quality substantially improves {xt =1-23). Except for the 
temperature (kT=36keV), some large differences are en- 
countered in the spectral parameters with respect to the 
BeppoSAX results. First of all the presence of the edge 
with Eo= 8.1 keV and an optical depth r — 0.13, which 
is not required, although not excluded by the BeppoSAX 
data. Here we note that although BeppoSAX has a bet- 
ter energy resolution than RXTE, the latter has a much 
larger collecting area. The edge energy corresponds, within 
the error, to Fe xix-xxii; the optical depth to an equiva- 
lent hydrogen column density of ^ 1.1 x lO^'^ApJ^ cm~^. 
The abundances results seem to be much lower than so- 
lar values. We note that fixing them to the solar value, 
the fit quality decreases noticeably (x^=2.13). The ori- 



gin of such large difference is unclear, but it should be 
noted that at the low resolution of RXTE the various 
iron lines are heavily blended. Also, the absorber is to- 
tal while column densities are within errors comparable 
to the BeppoSAX fits. The presence of such an absorber 
could be consistent with the presence of the absorption 
edge. The equivalent width of the K^ iron line (216 eV) is 
similar to that derived from the BeppoSAX data. With 
this model the observed fiux in the 2-10 keV range is 
6.64x10^^^ erg cm~^ s~^, ^ 1.4 times larger than that ob- 
served in July 1998 by BeppoSAX. An attempt to fit the 
RXTE data with an additional reflection component has 
been performed, resulting in a modest improvement in the 
quality of the fit, significant at the 95.9% confidence level 
(F-test; see Table 2 model 4 shown in Fig. 8). The temper- 
ature is lower (kT=26keV) than that derived in model 3 
and similar to that derived from BeppoSAX data. The 
relative normalization of the reflection component is com- 
parable, within errors, to the BeppoSAX fits. In summary 
the RXTE grand average spectrum provides the first evi- 
dence of absorption from ionized material. 

We also performed a phase-resolved spectral analysis of 
the RXTE data, limiting ourselves to study the maximum 
and minimum phases. Within the errors, no parameter 
changes between the two phases apart from the column 



density of the cold absorber, which is 1.37" 
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Fig. 7. Variations of the covering fraction of the partial 
absorber and E.W. of fluorescence iron line along the spin 
pulse period obtained from the BcppoSAX observation. 
A composite model consisting of a MEKAL emission plus 
interstellar hydrogen column density, partial absorber and 
Gaussian line has been fitted to the combined LEGS and 
MEGS spectra extracted in 10 spin bins. The temperature, 
the interstellar and partial absorber column densities and 
abundances have been kept constant (see text). 



5. Discussion 

The RXTE and BeppoSAX observations have allowed, for 
the first time, a simultaneous characterization of the X-ray 
temporal and spectral behaviour of V 709 Gas. 

5.1. Variability 

It has been possible to confirm the rotational period of 
the WD at 312.75s and to rule out the presence of any or- 
bital and sideband periodicities. This system is therefore 
a typical disc accretor. The rotational modulation is re- 
markably non-sinusoidal compared to most IPs. The spin 
pulses are clearly observed at all energies up to 25keV, 
with decreasing amplitudes from soft to hard X-rays, a 
typical behaviour of IPs. Furthermore, we have found that 
the power at the spin frequency and first and second har- 
monics changes with time: the spin pulsed fraction ob- 
served by BeppoSAX in July 1998 was twice that ob- 
served in February of the same year by ROSAT (compar- 
ison being done in the 0.1-2.4keV range) and about 1.2 
times larger than that observed by RXTE in March 1997 
( comparing the 2-5 keV range). Not only the amplitude 



but also the morphology of the rotational modulation has 
changed. The one observed by BeppoSAX is highly asym- 
metrical with a maximum possessing two peaks with dif- 
ferent intensities in the hard band, but it is single-peaked 
in the soft (ROSAT-like) band. The latter contrasts with 
the double-peaked maximum (with similar intensities) as 
observed by ROSAT. The BeppoSAX pulsation also dif- 
fers from the very symmetrical RXTE spin pulse which 
presents a relatively fiat~top maximum and more resem- 
bles the ROSAT light curve. It is worth noting that the 
pulse maximum extends over the same phase range in all 
observations. The presence of the first harmonic in the two 
observations in 1998, contrasts with the lack of detection 
in March 1997. A strong first harmonic, compared to the 
second, would produce a double peaked light curve with 
two maxima offsctted in phase by 180°. In V 709 Gas the 
second harmonic is the stronger and produces an asym- 
metrical profile when the first harmonic is present and 
strong, as observed by BeppoSAX. When the first har- 
monic is weak or even absent the curve possesses similar 
maxima as observed by ROSAT and RXTE. 
The changes in the shape of light curve are also accom- 
pained by variations in the X-ray fiux. Gomparing the 
ROSAT PSPG count rates observed in 1992 (Motch et al. 
1996) with those converted from the ROSAT HRI pointing 
in February 1998 (Norton et al. 1999), we find that V709 
Gas was brighter by a factor of ~ 1.2 in the latter obser- 
vation. Furthermore converting the BeppoSAX flux in the 
ROSAT energy pass-band, we find that the binary was 
brigher by a factor of ^1.7 in July of the same year. This 
suggests that V 709 Gas experiences changes in luminosity 
on time scales of less than months to years and supports 
the identification of a highly variable X-ray source made 
by Motch et al. (1996). Hence changes in the accretion 
rate are likely to occur in this system (see sect. 5.2) which 
give rise to changes in the amplitude and shape of spin 
pulses. This phenomenon has also been observed in other 
IPs such as FOAqr (Beardmore et al. 1998; de Martino 
et al. 1999), TXGol (Norton et al. 1997) and BGGMi (de 
Martino et al. 1995). As the accretion disc extends down to 
the magnetospheric radius, which scales as ex M^^/^, an 
increase in the mass accretion rate can produce a shrink- 
ing of the magnetospheric radius and hence an accretion 
spread over larger areas of the WD surface. At higher ac- 
cretion rates the effects of absorption are expected to be 
larger, thus influencing the shapes and amplitudes of the 
low energy rotational modulation. The increase in the spin 
amplitude between the February and July 1998 could be 
consistent with this interpretation. On the other hand, the 
March 1997 RXTE spin pulse has a lower amplitude with 
respect to that in 1998. At highest accretion rates, some 
of the local absorbers become ionized, which can reduce 
the low energy modulation amplitudes. 
The double-humped light curve observed by ROSAT was 
attributed to the contributions of both accreting poles by 
Norton et al. (1999), but with an offset of the dipole axis 
thus producing two maxima separated by less than half a 
spin period. In the two-pole accretion scenario there are 



de Martino et al.: The X-ray emission 



11 



two ways to produce a double-peaked spin pulsation. One 
invokes large areas, due to the low magnetic field of the ac- 
creting WD, with vertical optical depth (along field lines) 
lower than the horizontal optical depth. The other invokes 
both wide and tall shocks, but with the vertical optical 
depth larger than the horizontal, as in the classical ac- 
cretion scenario. The RXTE and BcppoSAX observations 
have shown that the spin pulse is energy dependent and 
the pulses are in phase at all energies (i.e. above and be- 
low lOkeV). To account for a rotational modulation above 
10 keV, reflection cannot be the main mechanism, since, 
as shown in sect. 4, it would produce an anti-phased pulse 
with respect to the low energy pulsation. A self occultation 
of tall shocks was proposed as a solution for the high en- 
ergy rotational modulation in IPs (Mukai 1999). This can 
account for the phasing of spin pulses at both high and 
low energies, for the lack of an energy dependence of the 
hard X-rays pulsations and for the sharpness of the rise 
and decay of the pulsation. Hence the possibility that in 
V 709 Cas the two accreting poles are wide and tall, with 
the optical depth along the magnetic field lines larger than 
across them (as in the classical accretion curtainmodel) , 
seems to be the most viable. Since V709 Cas is a classi- 
cal disc-accretor, the accretion is expected to occur onto 
both poles. Also, in order to have the same phasing at 
both high and low energies, i + S < 90°, where i is the bi- 
nary inclination and 6 is the magnetic colatitude (Mukai 
1999), with an upper limit to the inclination of z < 68° due 
to the lack of an X-ray eclipse in this system. This would 
imply that the upper pole dominates the hard X-rays at 
all phases and hence does not produce a strong modula- 
tion, but the lower pole produces a maximum when the 
upper pole points away from the observer. The lower pole 
appears into view for ~ 40% of a cycle, thus producing the 
sharp rise to and sharp decay from the maximum. When 
the upper pole points towards the observer (spin mini- 
mum), the lower pole is out of view. A simplified sketch 
of the geometry is shown in Fig. 9 for an illustrative set of 
values, i — 45°, S — 30° and shock height h = 0.2 i?^^ for 
two symmetrical poles. An asymmetry between the two 
poles was pointed out by Norton et al (1999) which seems 
to be necessary to explain the details of the light curves. 
An offset value of ^ 108° could be derived from the sepa- 
ration of the two maxima in the ROSAT HRI light curve. 
However for a detailed modeling of the accretion geometry 
one should take into account the whole three-dimensional 
shape of the emission region which would imply too many 
degrees of freedom. For our purposes, we only note that 
the net effects of the absorption within the two interven- 
ing accretion curtains at their respective spin maximum 
will then harden the spectrum at these phases. This ex- 
plains the lack of spectral changes during the maximum 
and dip phases. Furthermore, the BeppoSAX observation 
shows that this maximum is highly asymmetrical, which 
would suggest that the contribution from some parts of 
the accretion regions are not identical at some epochs. 
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Fig. 9. A schematic sketch of the X-ray spin modulation 
in V709 Cas for i = 45°, S = 30° and a shock height of 
0.2Rwd assuming for semplicity two symmetrical accreting 
poles. At spin minimum (the observer's line of sight is 
indicated by the arrow from the upper right), the lower 
pole is hidden, while the upper pole is seen along the field 
lines. At spin maximum (the line of sight is shown from 
the left) , the entire upper pole and most of the lower pole 
are visible (the regions hidden by the WD are indicated 
by the dashed lines), and both poles are seen across the 
field lines. The lower pole is visible for ~ 40% of a cycle, 
thus being responsible for the broad maximum. The upper 
pole does not produce a strong modulation, being visible 
at all spin phases. 

5.2. Spectral Properties 

The spectral analysis carried out from 0.1 to 100 keV has 
shown that V709 Cas is a classical IP in its X-ray char- 
acteristics. A model consisting of a single temperature 
optically thin emitting gas, attenuated by both interstel- 
lar material and a partial absorber, and an emission fea- 
ture corresponding to the fluorescent 6.4 keV iron line, 
describes well the BeppoSAX X-ray spectrum. The post- 
shock temperature ranges between 42 and 27 keV, depend- 
ing on whether reflection from the WD surface is included 
in the spectral fits. Within errors the difference between 
the two temperatures is not significant. We also note that a 
multi-temperature plasma is not required by either RXTE 
or BeppoSAX data. In this respect, we compare the de- 
rived temperature of the post-shock region with that ex- 
pected from the simple shock model 

kTshock = 3/8 G /^ mn M^d Rwd 
with fi the mean molecular weight and mn the mass of hy- 
drogen. From optical spectroscopy, Bonnet-Bidaud et al. 
(2001) constrain the WD radius at a value ^ 8.7 x 10^ cm 
and a mass of 0.75Mq, which then implies a temperature 
of 30keV, in close agreement with our BeppoSAX and 
RXTE results. 

From the spectral fits we have found that the column den- 
sity of the partial absorber cannot alone explain the large 
equivalent width (~ 200 eV) of the 6.4 keV line. Hence 
the contribution from reflection by the cold material of 
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the WD surface should dominate the fluorescent iron line. 
At a temperature of ^30keV, an EW > 100 eV is ex- 
pected for 6 < 36° (Matt 1999), where 6 is the viewing 
angle between the line of sight and the direction orthog- 
onal to the illuminated portion of the WD surface. From 
Fig. 9, the largest contribution to the iron line from the re- 
flected region of the WD surface should then be produced 
by the upper pole at spin minimum (9 = 15°). However 
if some asymmetry exists, this enables to see the lower 
pole more effectively, thus giving rise to a large equivalent 
width also at other spin phases. We have indeed found 
spectral variations with the rotational period, which indi- 
cate a spin variability of the partial absorber and emission 
line parameters, the largest contribution of both being at 
rotational minimum. Again, the accretion curtain scenario 
matches well the spectral behaviour of spin pulses. When 
the curtain points towards the observer, the absorption 
effects are larger (the covering factor is larger) as well as 
the effects produced by reflection, since the accretion area 
projected on the WD is larger. This is also seen from fit- 
ting the spectra at spin maximum and minimum which 
show an increase in the Compton reflection component at 
spin minimum. 

The spectral analysis of 1997 RXTE data has furthermore 
shown the presence of an absorption feature at '^ 8.1 keV 
which is interpreted as an absorption edge from Fe xix- 
XXII with an optical depth r ^ 0.11 — 0.13. This would im- 
ply an ionization parameter of ~ 100 (Kallman & Mc Cray 
1982). This is the first detection of an ionized absorber in 
an IP. Such material could be located in the pre shock 
regions of the accretion funnel. The lack of detection of a 
variability of this feature at the rotational period might 
suggest that this material is distributed over a large angle 
as seen by the observer. The results of the RXTE spectral 
fits also suggest the presence of a total absorber, which is 
consistent with a large covering area. 
We attempt to obtain estimates of the accretion rate 
from the X-ray luminosities of BeppoSAX and RXTE ob- 
servations. Bonnet-Bidaud et al (2001) constrain V709 
Cas at a distance between 210-250 pc, which would im- 
ply a bolometric luminosity of 1.1-1.5 xlO^'^ergs^^ dur- 
ing the March 1997 RXTE observations and of 0.8- 
l.lxlO^'^ergs^^ during the July 1998 BeppoSAX observa- 
tions. This gives an accretion rate of 9. 6-13. Ox 10^'"' gs~^ 
in 1997 and of 6.9-9.6xl0i5gs-i in 1998. The accretion 
rate has then decreased by a factor of ^^1.4 between 1997 
and 1998. Hence the detection of an ionized material in 
the RXTE data could be consistent with a higher accre- 
tion rate experienced in 1997. 

We also derive an average emission measure E.M.=2.6- 
3.8 xlO^^cm-^ and E.M.= 1.7-2.4 x lO'"^'"^ cm^^ for the 
March 1997 and July 1998 observations, assuming the 
above range of distances. As EM ~ nj^ x l^, we give an 
upper limit to the number density of hydrogen and an esti- 
mate to the linear dimension of the post-shock region. We 
use the RXTE results since from the iron absorption edge, 
we have derived the equivalent hydrogen column density 
Nh "^ 1.1 X 10^^ Ap^^ cm~^. Assuming that this material 



is in the pre-shock region, and for simplicity in the imme- 
diate vicinity of the shock, the density in the region just 
below the shock is a factor of four higher. Then, from the 
range of emission measure, we obtain / ^ 3.4 — 4.9 x 10^ cm 
and uh ^ 1.8 — 2.6 x 10^® cm~'^. The former gives a non- 
negligible shock height consistently with a tall shock as 
suggested from the spin light curve. 



6. Conclusion 

The BeppoSAX and RXTE observations have allowed for 
the first time the study of the broad band X-ray emis- 
sion spectrum and variability of V 709 Cas which can be 
characterized as follows: 

a) V 709 Cas is dominated by the rotational pulsation 
of the accreting WD at a 312.75s period, with no sign of 
orbital or sideband periodicities. This confirms that it is 
a disc accretor. 

b) Its spectrum is hard and well described by an 
isothermal optically thin plasma at 27keV with complex 
absorption and an iron Kq fluorescent line, due to reflec- 
tion from the WD surface. Evidence of cool and ionized 
absorbing material in the pre-shock region is found in this 
system, e.g. from a partial absorber and an iron absorption 
edge. 

c) The rotational pulsation is compatible with com- 
plex absorption dominating the low energy range, while 
at higher energies pulsation is likely due to the presence 
of a non-zero shock height above the accreting poles. 

d) Variations along the rotational pulse in the partial 
covering absorber and reflection are compatible with the 
classical curtain scenario, where accretion material flows 
from the disc towards the polar regions of the WD via 
arc-shaped curtains. 

e) The shape and amplitude of the spin light curve as 
well as the X-ray flux change with time indicating that 
V709 Cas experiences variations in the mass accretion 
rate on timescales of less than months to years. 
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